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^ _ We continue the study of the non-metric theory of gravity introduced in 

■ hep-th/0611182 and gr-qc/0703002 and obtain its general spherically symmetric vac- 

uum solution. It respects the analog of the Birkhoff theorem, i.e., the vacuum spheri- 
cally symmetric solution is necessarily static. As in general relativity, the spherically 
symmetric solution is seen to describe a black hole. The exterior geometry is es- 



o 



sentially the same as in the Schwarzschild case, with power-law corrections to the 

, Newtonian potential. The behavior inside the black-hole region is different from 

> 

' the Schwarzschild case in that the usual spacetime singularity gets replaced by a 

■ singular surface of a new type, where all basic fields of the theory remain finite but 

jy-^ . metric ceases to exist. The theory does not admit arbitrarily small black holes: for 

o , 

, small objects, the curvature on the would-be horizon is so strong that non-metric 



modifications prevent the horizon from being formed. 



^ . The theory allows for modifications of gravity of very interesting nature. We 

H ; 

, discuss three physical effects, namely, (i) correction to Newton's law in the neigh- 

borhood of the source, (ii) renormalization of effective gravitational and cosmological 
constants at large distances from the source, and (iii) additional redshift factor be- 
tween spatial regions of different curvature. The first two effects can be responsible, 
respectively, for the observed anomaly in the acceleration of the Pioneer spacecraft 
and for the alleged missing mass in spiral galaxies and other astrophysical objects. 
The third effect can be used to propose a non-cosmological explanation of high red- 
shifts of quasars and gamma-ray bursts. 
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I. INTRODUCTION 

This is the second paper in the series devoted to the detailed study of the non-metric 
theory of gravity introduced in As in the case of general relativity (GR), one of 

the first applications of any theory of gravity must be the description of the geometry of a 
strongly gravitating object — a black hole. This is the basic aim of the present paper, in 
which we obtain and analyze the spherically symmetric vacuum solution. The considerations 
of this work will also serve as an illustration to the rather abstract formalism of (2I]. 

The theory of gravity under investigation is defined by the action 

S[B, A, ^] = [ B' AF' + - + S'^) B' A B^ , (1) 

= [Tr (^2^ ,Tr (^^)] . 

Here, 5* is a (complex) su{2) Lie-algebra valued two-form (the indices i,j = 1,2,3 belong 
to the su{2) Lie algebra), = dA^ + |[y4,y4]* is the curvature of the 5u{2) Lie-algebra 
valued connection A*, is a traceless symmetric "Lagrange multiplier" field, and G is the 
Newton's constant. The function [Tr (\1'^) , Tr {^^y\ is a function of two scalar invariants 
that can be constructed from and is responsible for the departure of the theory under 
consideration from GR. The constant part of the function is (a multiple of) the usual 
cosmological constant A. 

The theory with action ([!]) modifies the Plebahski self-dual formulation of general rela- 

n , . 

tivity [3j, in which cp = (pQ = —A/3. In the Plebahski theory, the field B^, according to the 
equations of motion, can be canonically decomposed into products of tetrad basis one-forms, 
which then define a unique distinguished metric satisfying the Einstein equation. With the 
appearance of a nontrivial function 0(\1/) in ([1]), this property is no longer valid: the aris- 
ing metric is defined only up to a conformal factor, and it does not, in general, satisfy the 
vacuum Einstein equations. 

The nature of this modification of gravity is quite special and deserves a few words. In the 
numerous existing schemes of modified gravity, one usually deals with a metric theory and 
modifies the Hilbert-Einstein action by introducing extra degrees of freedom — either by 
increasing the number of derivatives (higher-derivative gravity), or by introducing extra fields 
(scalar- vector-tensor theories), or by considering extra dimensions (braneworlds) . Nothing 
of the listed takes place in our generalization: the theory remains four-dimensional, there 
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are no extra fields, the number of derivatives is not increased, and, as in GR, tliere are still 
just two propagating degrees of freedom, as can be seen from the canonical analysis of the 
theory; see, e.g., j4|. This is achieved by first recasting the theory in a form which does 
not contain any metric [3], and then modifying the theory in this non-metric form. It is 
thus important to stress that the term "non-metric" is used here not in the sense that some 
additional degrees of freedom are present along with the usual metric, but in the sense that 
metric does not even appear in the formulation of the theory. 

The difference between our modified gravity and a variety of other approaches can also be 
seen from the fact that our theory, while modifying the spherically symmetric solution, does 
not modify the formal general-relativistic cosmological equations. Indeed, due to the high 
symmetry, the field which is a close analog of the Weyl curvature spinor in our theory, 
is identically zero for cosmological solutions, which then coincide with those of general 
relativity with a cosmological constant determined by 0(0). However, any departures from 
homogeneity and isotropy will be essential, so that the evolution of perturbations in our 
theory will, in principle, be different from that in the concordance LCDM cosmology. This 
issue will be studied separately. 

A specific choice of the function (p uniquely fixes a theory from the class ([1]). One way of 
fixing the form of this function is to regard it as an effective quantum contribution to the 
original classical action. In this respect, it was argued in l| that the class of theories ([1]) is 
closed under the renormalization-group flow. The conjecture of asymptotic safety applied 
to the case at hand then asserts the existence of a non-trivial ultra-violet fixed point of 
the renormalization-group flow described by a certain function 0*(\I'). It would then make 
sense to choose this fixed-point function 0*(\I') in ([1]) because the quantum field theory so 
defined would have extremely appealing properties: the action would not get perturbatively 
renormalized, describing an essentially finite quantum theory of gravity. This gives a possible 
scenario for fixing from the theoretical side. 

In a pure quantum theory of gravity, the Planck scale is the only scale that can enter 
the renormalized action. However, once the theory is coupled to matter, the corresponding 
quantum loops in Feynman diagrams will also affect the function cj). Thus, it should be 
expected that the ultra-violet fixed-point function 0* will contain more than one physical 
scale. Qualitative features of the function 0* can be anticipated from the fact that, as we 
shall see in this paper, a theory with a non-constant function 0(\l/) in many respects be- 
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haves like a theory with physical parameters depending on the scale. Indeed, as we shall see, 
the function 0(\E') itself will receive an interpretation of a curvature-dependent cosmological 
"constant" . Another effect is that passing between spatial regions of different curvature gen- 
erally introduces a renormalization in the strength of the gravitational interaction. All this is 
strongly suggestive of the renormalization-group flow phenomena, where one has to identify 
\E' (having the dimension of the curvature) with the energy scale squared. The crucial differ- 
ence between our scheme and the renormalization-group flow familiar from the framework of 
effective field theory is that similar effects arise here in a diffeomorphism- invariant context. 
This renormalization-group interpretation suggests that the function 0*(\E') must look as a 
sequence of plateaus, with crossover regions between plateaus corresponding to the scales 
where new physics (new degrees of freedom) come into play. This discussion motivates some 
assumptions we make about the form of the function 0(\1') in the section where we discuss 
possible long distance modifications of gravity. 

Of special importance is the question of coupling our theory of gravity ([T]) to other fields. 



As was discussed in [2[, coupling to Yang-Mills fields (or electromagnetic field) is seamless, 
and the action described in [5] extends to the non-metric situation without any problem. 
This action tells us how massless particles such as photons interact with gravity and predicts 
their motion in a non- metric background. It is found that, in the approximation of geometric 
optics, photons move along null geodesies of the metric determined by the field B^. The 
ambiguity in the choice of the conformal factor discussed in detail in does not affect this 
conclusion because the paths of null geodesies (but, of course, not the affine parameter along 
them) are independent of this choice. An important issue that has not yet been addressed 
in the framework of non-metric gravity and that prevents us from a complete analysis of the 
physical predictions of theory ([1]) concerns coupling to (massive) matter degrees of freedom. 
The action for a massive field of spin 1/2 proposed in 5(] for the Plebahski formulation of 
general relativity turns out to be incompatible with the non-metric character of the theory 
under investigation, hence, calls for revision. 

Quantitative predictions of our theory, e.g., concerning the motion of stars and gas in 
galaxies, will depend on the details of the coupling of massive matter to the basic gravita- 
tional degrees of freedom. However, its certain general features can already be described in 
the absence of these details. For instance, in the domains of "metricity," in which the func- 
tion (f) is almost constant (where its dimensionless derivatives \d(j)/d'$\ <^ 1), the present 
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theory of gravity behaves very closely to general relativity, and it is quite reasonable to 
expect that matter will also behave accordingly, moving relativistically in the background 
of the arising metric. Assuming that several such domains of "metricity" exist at different 
scales of curvatures, we find that the effective gravitational mass of a central body is dif- 
ferent in the corresponding spatial regions. From this simple observation one can conclude 
that the effective gravitational mass of a body continuously depends on the distance to this 
body even in the case of general function — the effect of scale-dependence of the grav- 
itational coupling which was mentioned above. One can use this property to account for 
the phenomenon of missing mass observed in gravitating objects such as spiral and elliptical 
galaxies. Another interesting general prediction of the present theory is the appearance 
of an additional redshift factor between regions of different space-time curvature. As we 
point out in this paper, this effect can be used to account for the observed high redshifts of 
quasars and gamma-ray bursts. However, practical use of these features of the new theory 
to explain physical phenomena requires the specific knowledge of the underlying function 
together with the analysis of the physical content and interpretation of the theory. 
This will be the subject of the future work. 

The organization of the paper is as follows. To obtain a spherically symmetric solution, 
we first obtain an ansatz for the symmetric field B\ This is done in Sec. [Tll in which we 
also analyze consequences of the modified "metricity" equations. The field equations are 
obtained in Sec. Illli Solutions are analyzed and interpreted in Sec. IIVI Possible modifications 
of gravity are discussed in Sec. |Vl and possible physical effects in Sec. IVD Our results are 
summarized in Sec. I VI II In the appendix, we give a detailed proof of the static property 
of the spherically symmetric solution in the theory under consideration and discuss in more 
generality some details of the large- distance modifications of gravity. 

II. SPHERICALLY SYMMETRIC ANSATZ AND METRICITY EQUATIONS 

A. Spherically symmetric Lie-algebra-valued two-form 

The most general spherically symmetric su{2) Lie- algebra- valued two-form can be ob- 
tained from the condition that an SO (3) rotation corresponds to a gauge transformation. 
There exists standard technique in the literature allowing one to obtain the relevant expres- 



sion; see, e 



One gets: 



(2) 



B = J2i B't' = {(pidt Ad9-xi sin 9 dr A dcj)) 
+ (02 sin 6dt Ad(j) + X2dr A d6) 
+ (03 sin 9d6 Ad(j) + x^dt A dr) . 

Here, (t,r,9,(f)) is the standard set of spherical coordinates, and (pi, i = 1)2,3, are 
functions of t and r only. The symbols r* denote the su{2) generators r* = — (z/2)cr\ where 
cr* are the Pauli matrices. As it turns out, it is much more convenient to work not with the 
adjoint, but with the fundamental representation of su{2). This amounts to working in the 
spinor formalism. We used spinor formalism rather heavily in , and will continue to do so 
in this paper. To pass from the SO (3) form of the fields to their spinor representation one 
has to replace every lower-case Latin index z, j, . . . = 1, 2, 3 by a symmetric pair of unprimed 
spinor indices. Equivalently, every 5u{2) Lie-algebra-valued field gets replaced by a 2 x 2 
matrix-valued field. Thus, it will be convenient to rewrite the above ansatz for B in terms 
of the matrices 



(3) 

















(::)■ 










which are related to r* as 


ri = 


(l/2z)(X+ 






have: 



B = X_B_+ X+B+ + XBo 



= X_ 

+ x+ 



^Y^cpidt + ]^X2dr^ AdO + {^^Xidr + ^02^?^ ) sin 61 A 
^ \dt + \x2dr \ AdO + ( -^Xidr + l-(f)2dt 1 sin 6^ A 



(4) 



2r 2' J V 2r 

- iX (03 sin 9d9 Ad(t) + x^dt A dr) . 
It is not hard to show, and this fact was used heavily in that, by a convenient choice 
of the spinor basis, the spinor counterpart of the quantity \E' can be put into the form 

= a[x^®X^ + X+® X+) + /3 ® X_ + X_ O X+ + 4X (g) . (5) 

The functions a and (3 are related to the invariant characteristics of ^ : 



Tr (^2) = 2a^ + 6/5^ , Tr (^^j = 6/3 (a^ - /J^) 



(6) 
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In view of ([6]), one can regard [Tr , Tr as a function of a and /3. 
In the spherically symmetric case, the field has the form 

= ^(r) (^x'x' - ^^''r') , (7) 

where the natural Euclidean coordinates realizing the group of rotations, 

and = (x*)^. This condition implies a = in ([5]), and the field ^ is parameterized by 
a single function /5 : 

= j3 (^X+ ® X_ + X_ ® X+ + AX ® . (8) 

In other words, the Lagrange multiplier field must be algebraically special, of type D. This 
is, of course, exactly the property of \1/ in a spherically symmetric solution of the usual GR. 
This property remains unchanged in the non-metric theory of gravity under consideration. 



B. "Metricity" equations 

With the above ansatz for the B field, it is easy to compute the quantity B^ A B^ that 
appears in the metricity equations — the equations stemming from ([T]) as the field \1/ is 
varied. These equations are discussed in [2] in great length, and we will not repeat that 
discussion here. We just note that, as it is easy to check, the quantity 5* A B\ with B^ 
given by ([2]) above, is diagonal as a 3 x 3 matrix. This implies that the matrix $ := 90/9^, 
which is the traceless part of 5' A B^ in view of the metricity equation, is also diagonal. 
Hence, 0q, := d(l)/da must be identically zero in view of equation (26) of [2]. The property 
a = ensures this condition, in particular, if the function is a regular function of its 
arguments ^ in the neighborhood of zero. 

Let us now write the metricity equations specialized to the type D at hand. They are 
easy to obtain from equations (29) of j2j| by setting 0q, = 0. We have 

5+ A 5+ = 5_ A 5_ = 5+ A 5o = 5_ A 5o = , 

2B+ A S_ + Bo A Bo = -20/3 [b+ A B_ - ^Bq A Bq 

where := d(j)/d(3. Let us now see what this implies about our ansatz (jl]). The last two 
equations in the first line of ([9]) are automatically satisfied, while the first two equations 
imply 0iXi = 02X2- This allows us to write the X± components in (jlj) as 

B_ = cm Al , B+ = cn Am , (10) 



8 



where 

V 02 / V 2? 2 y ^^^^ 

:= ^ (^-^d9 - i sine dcp^ 



m := ^ [ —^d6 — isinO d(f) ] , n := rj ^——(l)2dt — —Xidr 

In these expressions, ^ and 1] are arbitrary functions of t and r, and c = ^rj. To fix these 
functions, we equate the component Bq in (jlj) to the third canonical two-form: 

- i (03 sin 9d9 Ad(f) + x^dt Adr) = lAn-mAm, (12) 

which gives and t]"^ in terms of 0j and Xi • 

201 02Xl 

Eventually, the expression for B takes the form 

B = X_c mAl + X+c nAm + X{lAn-mATh). (14) 

The last metricity equation [the second line of (Q] then relates the function c to the function 
/3 in ([8]) through a derivative of : 



c^^l^. (15) 

Now we are going to simplify the expressions for one- forms (ITT]) . First, we can choose ^ 
as a new radial coordinate. After this, introducing new functions /, g, and h, one can write 
the one-forms /, n, m, and rh as 

1 1 _ r 

I = --j= {fdt — gdr) , n = — ^ {fdt + gdr) , m, m = — ^ {hd9 ± z sin ^ dcj)) . (16) 

By solving the system of field equations, one can prove that the function h is just a 
constant, and that the sought functions /, g, and (3 are independent of time. We demonstrate 
this property in the appendix. Then, by rescaling the angle and the radial coordinate r, 
we can set h to be identically equal to unity, after which the canonical set of one-forms is 
expressed as 

1 1 _ r 

I = -j= (fdt — gdr) , n = — ^ (fdt + gdr) , m, m = — ^ {d9 ± i sin 9 dcj)) , (17) 

and the metric ds^ = 21 ®n — 2m ® fn defined by tetrad (ITTI) assumes the standard form 

ds^ = fdf - gHr^ - r^ {d9' + sin^ 9 dcj)^) (18) 
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of the spherically symmetric problem of general relativity. In this form, we assume all 
space-time coordinates to be real. Thus, the only novelty as compared to the GR case is 
the presence of the function c in the first two terms of ( JT4l) . The function c is related to 
P through the derivative of via ( |T5l) . Note that, if 0^ = 0, we have c = 1, and 

the two-form B reduces to that of the GR (with a cosmological constant determined by the 
constant value of 0). 



III. FIELD EQUATIONS 

The theory under consideration respects the analog of Birkhoff's theorem. The proof 
of the static property of the metric is given in the appendix; its asymptotic flatness is 
demonstrated below. In the main text, we just assume the tetrad forms to be given by 
expression (ITTj) . and the field by expression (IHl), with f,g, and (3 being functions of the 
radial coordinate r only, and we will be looking for solutions of these functions. 



A. Structural equations 



The Cartan structural equations (called compatibility equations in [2]) give an algebraic 
relation between the connection form A and the two-form B and its exterior derivative, 
which allows one to solve these equations with respect to A. As the first step towards this 
solution, we obtain manageable expressions for the quantities dBj- and dB^, where 



B_ = cmAl, B^ = cnAm, Bq = I A n — m A m . 



(19) 



and the one-forms /, n, m, m are given by (fT7|) . This is an easy exercise in differentiation, 
similar to what one does in obtaining the Ricci rotation coefficients for metric f|T8l) . After 
simple calculation, we obtain 



dB^ 



dB, 



(rc/)' 



gf 
gf 



I An Am — c cot 91 Am Am 
I A n A m + c cot 9 n Am Affi 



(20) 



V2 



dBo = (I — n) Am Am . 

rg 



where the prime denotes differentiation with respect to r. 
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Another way of expressing the results for dB^, dB^ is to project the arising three-forms 
onto the basis of dual one-forms. This is conveniently done by introducing, for an arbitrary 
four-form C, the scalar quantity defined as 

C = C^l An Am Am. (21) 

We get 

(rcf)' 

{dB_ A m)^ = {dB+ A m)^ - ^ ^ ^ 



V2rgf' 

- {dB_ A n)^ = {dB+ A l)^ = cot 9 , (22) 

a/2 

{dBo A n)^ = {dBo A l)^ = — , 

rg 

with all other components being zero. 

The Cartan structural equations that determine the components A±, Aq of the connection 
have the form 

C^ = A_ ABq- AqAB_, 

C+ = AoA5+-A+Afio, (23) 
Co = A_ A fi+ - A , 

with 

C± = -dB±, C, = -Ub,- (24) 

see {2! for derivation. These equations are obtained from the equations in by specializing 
to the case of a "constant" basis in the space of "internal" spinors, see for a description of 
the distinction between "internal" and "spacetime" spinors. The basis in which the Lagrange 
multiplier has the form ([8]) is precisely such a basis. 

These equations are solved by first computing their components using definition (pTI) . 
which gives 

(C_ A /), = (A_)„ , (C_ A n), = -{A_)i + c(Ao)^ , 

(C_ A m)„ = -{A_)fn , (C_ A m)„ = {A_)^ - c(Ao)„ , 



(C+ A l)^ = c{Ao)m - {A+)n , (C+ A n)^ = {A+)i , 

(C+ A m), = -c{Ao)i + {A+)rn , {C+ A m), = -(A+)^ 

(Co A/)^ = c(A_)„, {CoAn)^ = c{A+)rn, 

(Co A m)„ = -c{A^)i , (Co A m)^ = -c(A+), 



(25) 
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and then solving this hnear system of equations, with the result 



(A„)„ = c-i(Co A/),, 

{A+)i = (C+ A n), , 
{A+)m = Am)^, 

{Ao)i = c-\Co A n)^ - c-\C+ A m)^ 



(A_)„ = (C_A/),, 
(A_)™ = Am),, 

(A+)„ = -c-\CoAm)^, 
{A+)fn = c-\CoAn)y, 

(Ao)„ = c-2(Co A - c-^(C_ A m). 



(26) 



(Ao),r^ = -c~^(Co A m)„ + c-i(C+ A /), , (Ao)^ = -0-2(^0 a m), + c-i(C_ A n)^ . 

Here, the subscripts /, n, m, and m indicate the corresponding components in the develop- 
ment over the basis one- forms (fTTIl . 
Using (I22I) and ([21]), we finally have 

1 . 1 



^0 



\plrcg 
{rcfy 



m . 



A, 



V2 



rcg 



■ m , 



V2 



9 L 



rcf 



1 



rc^ 



ri \ coiO _ 

[I + n) ■;=— [m — m) 

V2r 



(27) 



This is our final expression for the spin coefficients. The usual metric case is obtained by 
setting c = 1. 



B. Field equations 

Now we are in a position to derive and solve the field equations of our theory which are 
analogs of (vacuum) Einstein equations in GR. The vacuum field equations are obtained by 
varying action ([T]) with respect to B: 

F+ ($ + 0Id)5 = 0. (28) 

The components of the curvature of connection (127|) are given by 

F_ = dA_ + A„ A v4o , F+ = dA+ + A , Fq = dAo + 2A_ A A+ . (29) 

Using the explicit spinor form ([8]) for we get 

(^ + 0Id)S = X.{p + 0)fi_ + X+(/3 + 0)S+ + X(0 - 2p)Bo . (30) 
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Thus, we can write ( 128|) in components: 

dA_ + A^AAo + {P + = , 

dA+ + AoAA+ + {p + 0)5+ = , (31) 
dAo + 2A_ A + (0 - 2p)Bo = . 



All computations are straightforward in view of (jlQll and (l27l) . Each of the first two equations 
in ( pTl) gives rise to the two differential equations 

^* ' ^l-ili = -^ = /5 + '^' (32) 



while the last equation in fl^ gives two additional equations 

1 



2/3, (33) 



rg* c2 ' 



where 

/*:=c/, g*--=cg. (34) 

In view of (1151) and by virtue of the Bianchi identity VF = 0, where V is the A-covariant 
derivative, only three of the four equations fl32|) and fl33l) are independent (this can be verified 
directly), and the system of three independent equations can be presented in the following 
convenient form: 

(2-0.)/?' = -^, ,- = l-(2/3-0)r^ ^ = ^^. (35) 
r f^g^ r{2 - 0/3) 

IV. SOLUTION AND ANALYSIS 

In this section, we proceed to solving the main system of equations f[T^ . fl5^ . which 
completely determine the two- form B given by f|T4l) and the one-form A given by fl27j) . 

A. Corrections to the metric case 

If our solution happens to be in the regime where |0/3| <^ 1, we have ~ 1, which is 
thus an approximate "metricity" regime. Equations (l35l) in this approximation lead to the 
Schwarzschild-(anti)-de Sitter form for the functions / and g in metric (ITSl) : 

/^ = ^, f = ,- = l-^ + 0or% (36) 
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where 0o is the value of the almost constant function (f){(3) in this domain, and is the 
Schwarzschild radius, which appears as the integration constant of the solution of the first 
equation in (l35l) . The constant 0o corresponds to the effective cosmological constant: 0o = 
-A/3. 

The approximate solution (1361) will, in particular, hold in the asymptotic region of large 
r if is an analytic function of its arguments Tr and Tr in the neighborhood of 
zero. In this case, as can be seen from ([6]), 0(/3) admits an expansion in powers of and 
f3^, and the leading term in this expansion for small (3 is 

m = <l^o±i'/3' + 0{(3') , (37) 

where £ is a constant of dimension length. Thus, we have 

0^ ^ ±2fp ^ ±^ (38) 

so that 

10^1 <1 for r>(£V,)'^^ (39) 

By solving fl5Bl) perturbatively in the small parameter i'^Vg/r^, one can obtain the next 
correction to solution (15^ in the domain of large r : 

"^^sO^S?)- 9" = l-7(l±Sf)+^«''- /;=C^(l=F^). (40) 



According to (fT5l) . the value of the "nonmetricity" parameter in this approximation is 
given by 

Thus, we can see that, under the assumption of regularity of the function 
[Tr (\E'^) , Tr (\I^^)] in the neighborhood of \1/ = 0, the solution in the asymptotic region 
r ^ oo tends to the metric form (the "nonmetricity" parameter c rapidly tends to zero) 
asymptotically describing a space of constant curvature. It is in this sense that an analog of 
the Birkhoff theorem holds in the theory under consideration. 

If £ ^ Ts, then the approximate regime fHOj) and fj4T|) is valid up to the "horizon" r = Vh, 
determined by the condition = g^"^ = 0. At the "horizon," the function c remains finite, 
and the function g^,, hence, also g, diverges. As in the metric case, this can be regarded as a 
coordinate singularity which can be removed by choosing new time and radial coordinates. 



14 



In this way, one can pass to the "black hole" region r < r^, in which the functions and 
g^^ change sign and become negative. Somewhere in this region, the condition |(/)^| <^ 1 
(or, equivalently, i'^Vg/r^ <^ 1) may cease to be valid, and the solution becomes strongly 
non-metric. To see what can happen in this region, we consider general solution. 



B. General solution 



First of all, one can note that the first and third equations in ( l35ll are singular at the 
point where 0/? = 2. One can remove this singularity by passing to a new radial coordinate. 
The value of (3 itself can be chosen as such a coordinate. Doing this, one can rewrite the 
system of equations (135|1 in terms of this new coordinate: 
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'2(3 



(42) 



This system is nonsingular and can easily be integrated. The function c^, as usually, is 
given by f|T5l) . Note that (3 has a physical meaning being a scalar characterizing the field 
according to ([6]) and (IHl). 

The metric flTSl) in the new coordinates {t, (3) is written as 



ds^ 



1 - 



-1 



fide - 1 



m 

3/3 



gld(3' 



\13) {de^ + sin'' e dcj)^) 



(43) 



where we have used f[T5|) and the first equation of fl42|) . 

Metric in this form indicates that the points where 0/? = 2, corresponding to = 0, are 
hypersurfaces across which the coordinates t and (3 change their space-time roles. At these 
hypersurfaces, metric fj43|) is degenerate. One can see that our basic fields remain finite at 
such points. Indeed, in the new coordinates {t,(3), the forms B± are proportional to the 
functions cf = and eg = g^,, which remain finite and nonzero at these points, and the 
potentially dangerous / A n component in Bq is also finite: 

dv T 
lAn = fgdt A dr = —c-'f,g^dt A = -— (1 + ^fs) f^gjt A d(3 , (44) 

where again we have used (fTSjl and the first equation of (142|) . The same is true for the spin 
coefficients A± given by f l271) . which are proportional to g~^. The potentially dangerous part 
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of the spin coefficient Aq is finite as well: 



i^ + ^) = ^(^-f + ^-^]dt = rf.g. (/? + (P)dt , (45) 



jrcf)' 1 
y/2g [ rcf rc^\ ' ' rg^. \ f* 

where we have used equation (!32|) . Due to Eq. (!28|l . the components of the curvature F are 
also finite. 

At a point where 0^ = —1, the quantity turns to zero. However, by similar reasoning, 
one can see that all components of A and B remain finite in the coordinates (t, (3). Therefore, 
the general solution is nonsingular at this point as well. The metric in the form fH3|) is 
degenerate, but this is not surprising as we are dealing with an intrinsically nonmetric 
theory, and the "nonmetricity" parameter c at this point is infinite. 

As we already noted in the previous subsection, the condition g'"^ = 0, hence = 0, 
corresponds to a horizon in metric (H3l) . This is a coordinate singularity which can be 
removed by passing to a Kruskal-like coordinate system. 

For an illustration, let us consider the function (f){(3) exactly in the form 

m = <l>o±i'P\ (46) 
with two possible signs. In this case, the solution is 

Af3) = ^ e±^^^ , g:' = 1-{2(3t fp' - <Po) r\f3) , f.g. = e^^'^ , (47) 



and the function c is given by 



= '^'^ (48) 

For i'^jS -C 1, this corresponds to the approximate solution (l40l) . (jH]). 

With the upper sign in (146|) . from the first equation of fHTl) . one can see that there is a 
minimum value of the radial coordinate 

rl = li'r,, (49) 

which is precisely the point where 0/3 = 2. 

In the case of the lower sign in (H^ . the function f3{r) is monotonic, and /5 — oo as r — 0, 
which resembles the behavior inside a classical black hole in GR. The point = 2 is absent 
in this case, but we have a singularity in metric (H3l) at the point £"^(3 = 1/2 corresponding 
to the condition (pp = —1. As we noted above, the fundamental fields A and B are regular 
at this point. 
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C. Minimal black hole 

For both signs in (H6ll . there is a minimum value of the "Schwarzschild radius" = 
for which the black-hole horizon exists. The mechanism, however, is different for the two 
signs. For the upper sign (corresponding to positive 0), there is a minimum value (H9ll that 
the coordinate r can take. This value of r is also the minimum of the function (yf~^. For 
the horizon to exist, this minimum value should be negative. Neglecting, for simplicity, the 
value of we can translate it into the condition 

rs>-J (50) 

for the Schwarzschild horizon to exist. If this condition is violated, then, instead of the 
Schwarzschild-like horizon, one has a naked surface of "non-metricity," which we describe in 
the following subsection. 

For the lower, negative sign in ( H6l) . the function g~'^ similarly has a minimum at /3 = 1/^^. 
The condition that this minimum is smaller than zero translates into the condition 

for the existence of horizon. 



D. Conformal structure 

In the modified theory of pure gravity described by action ([T]), the notion of a unique 
distinguished metric is replaced by the conformal class of metrics with respect to which the 
two-form B is self-dual, with metric ( IT8|) being a representative of this class. A distinguished 
physical metric from this class may arise only after one considers gravitational interaction 
of matter and radiation. Several potential candidates for such a metric can be envisaged at 
this level. Given the spinor two-form B, one can consider the Urbantke metric (7^^, defined 
by the relation [3] 

VW\9lu = ^e^^^^Tr {B,^Bp,Bs.) , (52) 

where the trace is taken with respect to the spinor indices. It is easy to see that the Urbantke 
metric is related to the metric defined by ( JTSl) by the conformal factor c^/^ . 
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Another distinguished conformal factor for the metric is obtained by the requirement that 
the quantity 



coincide with the volume element defined by the new metric. The metric line element that 
arises this way is given by 



where ds^ is the line element given by ffTSjl . 

In the metric case, where c = 1, all such definitions coincide with our "canonical" metric 



As we noted, before one considers coupling of our theory of gravity to matter, it is 
impossible to distinguish any of the listed possibilities for the metric as being the physical 
one. However, a rather natural requirement that the coupling of matter degrees of freedom 
to B is at most quadratic in B (the coupling of YM fields satisfies this requirement) favors 
the metric defined by (15^ . Indeed, the usual matter coupling to the Urbantke metric ( 152|1 
would be non-polynomial in the B field, which is undesirable for many reasons. 

In spite of the mentioned ambiguities, the presence of a distinguished conformal class of 
metrics, with respect to which the two-form field B is self-dual, allows us to speak about 
the conformal structure of the obtained solution. This conformal structure can be shown to 
have physical meaning reflecting the geometry of propagation of light. 

The conformal structure of our black-hole solution in the (t, (3) coordinate plane will 
strongly depend on the shape of our basic function 0(/?). However, its key details are easy 
to understand by looking at the form of metric fj43l) in the (t, (3) coordinates. The metric 
described by (143|) has the following types of critical surfaces defined by special positions in 
the [3 coordinate: 

1. The point where g'"^ = 0, hence, also = [the product f^:g^ is positive and finite 
for finite /? in view of the last equation in fH2l) ]. It defines a null horizon, analogous 
to the Schwarzschild horizon, separating different space-time regions in the black-hole 
solution. This is a coordinate singularity in (t, r) or {t, (3) coordinates, which can be 
removed by proceeding to Kruskal-like coordinates. 



-Tr(i?A5) 



(53) 




(54) 
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2. The point where 0^ = 2, or (? = 0. By differentiating the function g~'^, one can verify 
that this is a critical point (maximum or minimum) of this function, as well as of the 
radius r(/3) as a function of (3. Hence, typically, this critical point will not coincide 
with the previous one, where g'"^ vanishes. This is a position of true singularity in 
the class of metrics (H3l) . However, as we said before, the solution in our basic fields 
A and B can well be extended beyond this surface. What happens at this point in 
terms of our basic two-form field B is that the components B± and Bq no longer span 
a subspace in the space of two-forms which is self-dual in any metric. One can notice, 
for example, that cgdr = g^ {dr/dj3) dj3 = ai this point, so the one-forms 

cn = d = --j=f^,dt , (55) 

and the anti-self-dual two-forms 

B- = cmAn, B^ = clAm, (56) 

coincide with their self-dual counterparts B_ and B^, respectively. This means pre- 
cisely that there is no metric with respect to which the tripple B±, B is self-dual. 

3. The point where (j)/3 = —1, or c^^ = 0. This is another singularity in metric (H3l) . At 
this point, we have / A ri = in view of ( l44l) . and the two-form Bq coincides with its 
anti-self-dual counterpart Bq : 

Bq = I /\ n — m /\ fh = — {I /\ n + m /\ fh) = Bq . (57) 

Again, the fields A and B are well-behaved at this point, so we can cross it and proceed 
to a neighboring space-time region. 

As an example, in Fig.[T]we have pictured the conformal diagram obtained for the function 
0(/?) = which has only critical surfaces of type[T]and[2l The Schwarzschild-like horizons 
are pictured by solid lines, while the dashed lines correspond to horizons of type[2l The radius 
r has an absolute minimum value, in this case given by (H9l) . which is reached precisely at 
horizons of typeO As /? ^ oo, which corresponds to r ^ oo, one approaches a "singularity" 
indicated by thick dashed lines. 

Since, in the absence of matter couplings, the physical metric is not specified, the issue 
of "geodesic completeness" of solution is not well defined in the purely gravitational theory 
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FIG. 1: Conformal diagram of the spherically symmetric solution with the function (/>(/?) = £^/3^ 
and with the Schwarzschild radius satisfying (left image) and violating (right image) condition 
()50p . Different regions are numbered in such a way that the coordinate t is timelike in the odd 
regions, and spacelike in the even ones. Solid lines between regions indicate Schwarzschild-like 
horizons of typedl at which g~'^ = 0. Thin dashed lines indicate singular surfaces of type [21 where 
metric ceases to exist. The flow of time is vertical in the white regions and changes to horizontal 
in the grey regions. Thick dashed lines indicate the true singularity, where r = and (3 = oo. The 
configuration on the left image extends periodically and indefinitely upward and downward. We 
are living in one of the regions of type I; the asymptotic spatial infinity in this region is denoted 
by , and the future and past null infinities are denoted by and , respectively. 

under consideration. Nevertheless, it is interesting and instructive to see whether the metric 
defined by (fT8|) is geodesically complete as one approaches the singularity /3 — oo in the 
regions of type VI and VIII. It is easy to see that null geodesic completeness in such regions 
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is equivalent to the divergence of the integral 



oo 



/ 



1 + 0/3 
/3 



f.g.r{P)dp. 



(58) 



Thus, solution (l47|l . fj48|) with the upper sign has convergent integral (158|) . hence, the metric 
( IT8l) defined by this solution is not null geodesically complete. However, it is timelike 
geodesically complete. Indeed, in the regions in question, /? plays the role of the time 
coordinate since ^ — 1/2 as /3 ^ oo. To the leading order g'"^ ~ P(3'^r'^{(3) and dr ~ 
{i"^ /3)r{P)dp. This gives g'^dr'^ ~ —2(?dp^/9P'^, and the proper time to infinity /5 — >• cxd is 
logarithmically divergent. 

For the lower sign in fj471) . fj48|) . integral fl58|) is divergent; therefore, the corresponding 
metric ( JTSl) is null geodesically complete in the region (3 — > oo, or r — > 0. However, it is not 
timelike geodesically complete since the timelike distance to the singularity J gdr oc j dr is 
finite (the coordinate r is timelike and r ^ in this case). 

The described behavior makes the "singularity" at /5 — oo rather an interesting place. 
We have depicted it as a null surface because, in some respects, it is reminiscent of the usual 
null infinity. The reader should, however, keep in mind that the detailed structure of this 
"singularity" is quite unusual, and strongly depends on the details of the behavior of the 
function (p a.s (3 oo. In contrast, the structure of the "non- metric horizons" at 0^ = 2, —1 
is universal. It is also worth emphasizing that the singularity at — > oo is located "outside" 
of the black hole on the left image in Fig. [1] in the sense that another Schwarzschild-like 
horizon has to be crossed to reach it. It is in this sense that the theory under consideration 
"resolves" the singularity inside the black hole. 

In the examples considered above, we have assumed that the function 0(/3) behaves as 
±£^/3^ as /3 — s> oo. Some other choices are possible and may be interesting to consider. Thus, 
it is easy to devise 0(/3) so that there is a maximal possible "curvature" in the theory. To 
achieve this, one can choose 0(/?) diverging for some finite value of /3. A good example of 
such a function is given by 0(/3) = a;^/£^(l — x^), where, as before, x = £"^(3. It diverges 
at the value Pmax = 1/^^, which plays the role of maximal curvature in the theory. The 
qualitative behavior of the conformal diagram for the spherically symmetric solution will be 
the same in this case. The details of the behavior at the singularity at jSmax will, however, be 
different. Another interesting possibility not considered in this paper is that, for Planckian 
curvatures, the function 0(/3) changes rapidly — so that, in particular, one of the non-metric 
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horizons at 0/3 = 2, —1 is reached — and approaches a large constant value for (3 much larger 
than Planckian. In this scenario, we would get a conformal diagram consisting of elements 
similar to the ones described above, with a metric universe with large cosmological constant 
on the other side of the black hole. It could be interesting to study all such possibilities in 
more detail. 



V. MODIFICATION OF GRAVITY AT DIFFERENT CURVATURES 

The theory under consideration leads to a possibility of scale-dependent modification of 
gravity of very interesting nature. 

Defining the effective Schwarzschild radius rs{P) by the relation 

r,{P) = 2Pr\P), (59) 

we can present our main system of equations fj42l) in the form 

d(3 2/5 ' ^* ^ ^ + ' d(3 2(3' ^ ^ 

Integrating the first and the last equations in ( l60l) . we get the following relations, valid for 
arbitrary Pi and P2 '■ 

rsiPi) f*{f3i) g*{(32) 

where 

Zi(3i,l32)=exp( [ ' ^dp] . (62) 



The first equation in fl6T|) shows that we are dealing with a theory in which the effective 
Schwarzschild radius becomes distance-dependent, and the second equation indicates the 
presence of an additional redshift /blueshift factor in the metric. 

Assume that there exist two domains in the space of values of /3, in both of which one 
has 10^1 <^ 1, hence, both of which are characterized by the condition of metricity. Then the 
Schwarzschild radii will be constant in the corresponding regions of the radial coordinate 
and will be related by (!6T|) in which /?i and /?2 are the representative values of /? in the 
corresponding domains. The time-time component of the metric, according to the second 
relation of (16T!) . will exhibit the relative redshift factor between these two regions 

(2) (1) 
900 _ Qrr „_2 

— - . [b6) 

900 9rr 
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The cosmological constant, in general, will also be renormalized. 

As an illustration for these effects, we consider an explicit example, which is defined by 
the function 

0(/3) = 0o-0ilog(l + x2) , x:=e^(3. (64) 
The derivative of is given by the expression 

2ax 



where 



a 



(65) 



(66) 



is a dimensionless parameter. Function fl64j) has the property that |0/3| <C 1 in the two 
asymptotic regions x ^ (1 + Ittl)""*^ and x ^ 1 + so that the theory is approximately 
"metric" in both these regions. 

Equations can be integrated in this case, with the solution 



2r3 



„ _« arctan x—^ 



f asuctanx—^ 

j*y* — fc^ 



(67) 



We have chosen the integration constants in fl67j) so as to obtain the standard 
Schwarzschild solution at small radial distances, ^ (1 + \a\)^^ i'^r^, where x ^ (1 + |a|). 
At these distances, we can expand the right-hand sides of solution fl67|) to obtain 

P2^ 



X 



1-- + (ax~') 
X ^ 



2r3 

We also have 



f^g^ = l-- + 0(ax-^) , x>(l + |a|) 

X 



3a 



= IH hC(ax~^), x>(l + |a|) 

X 



(69) 



Collecting all the terms, we obtain the leading contribution to the Qqo coefficient in met- 
ric (fTSi) and ^00 coefficient in the Urbantke and "volume" metric (which coincide in this 
approximation) : 



5^00 



fi'oo 



1- - + 

r 



1-^ + 
r 



- 201 1 + log 



-201 



log 



2r3 



r — 



£V7 



(70) 



1 + \a\ 



At large radial distances ^ {1 + \a\)i'^rs, where x <^ (1 + |a|) , we have 



2r3 



e ''■2' X [l + ax + O (ax^)] , f^g^ = e [l + ax + O (ax^)] 



(71) 
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= l + 3ax + (ax^) , x < — \— . (72) 
^ ^ 1 + |a| 

In this case, we obtain the following leading contribution to the corresponding metrics, up 
to terms of order i'^Vs/r^: 



900 = e-"" 



^00 = e-™ 



1-62— + 4>or^ - e 2 — - 
r 2r 



7Va V a 



1 — 62 — + (p^r"^ 



(73) 

> (1 + \a\)l^rs 



s ■ 



r 

We note that ~ 1 in our approximation in these regions, so the physical metric is well 
defined. In fact, the contributions proportional to a in (170|) and (173!) are small and can be 
dropped if a is not very large. 

In the intermediate region of radial distances 

1 ^ 



< 



<l + l«l, (74) 



1+ |a| ~ 

the metric coefficient behaves in a complicated way; moreover, the physical metric is not 
well defined in this region unless a <^ 1. However, even without the detailed knowledge of 
matter couplings to gravity, simply on the basis of physical continuity, it is obvious that 
the apparent mass (or gravitational coupling) will be varying continuously between the two 
asymptotic regions. The intermediate region fTMl) is rather extended if \a\ ^ 1. 

It is interesting to note that the shift in the cosmological constant between the two 
asymptotic regions will be absent in a special case where the function 0(/3) takes the same 
value in the corresponding domains of the (3 space. In spite of the absence of such a shift, the 
Schwarzschild radius will still get renormalized. This can be illustrated by another example: 

2 

(/.(/?) = (/.o- (/.I x:=ep, (75) 

which is a function of both invariants Tr (\E'^) = 6/?^ and Tr (\E'^) = —6/5^ in the case a = 
[see Eq. (|6])]. This function is characterized by the property that (p ^ (pQ in both asymptotic 
domains x <^ 1 and x ^ 1. The mass renormalization exponent in (]B9P is given in this case 
by the integral 

oo ^ 



where a = I'^c 



'1- 



24 



In the case under consideration, we obtain the following approximate expressions for the 
metric, respectively, at small and large radial distances: 



r., , 2 ~ / 2r 



fl-oo = 1 ^4>or +2a 

r \ I'^r. 



3 \ 2 



c2 1 '"^ I J, 2 3 / 2r 



3 \ 2 



1 + \a\ 



(77) 



2Tra. 

goo = e 



1 — e3v^ h 00^^^ 

r 



e3v/3 - 



2^ 



27ra 

^00 = 6 3^ 



1 — 63^ \- (j)^r 

r 



> (1 + \a\)frs . 



(78) 



More general set of examples, characterized by arbitrary powers in the region of large and 
small curvatures is considered in Appendix [Bl The above formulas are in correspondence 
with the generic expressions flB4p and fIBSp . 

VI. PHYSICAL EFFECTS 
A. Modified gravity instead of dark matter 



Modifications of gravity at large distances are currently under consideration as an alter- 
native to the dark-matter phenomenon (see, e.g., Sj). Our theory can be a viable candidate 
in this respect. Indeed, we have seen in the previous section that the gravitational strength 
of the central spherically symmetric body is distance-dependent. Consider, for definiteness, 
our example Because the parameter a stands in the exponent of the renormalized 

gravitational mass [see Eq. fl73l) ]. one should have positive a ~ 1 in order that the gravita- 
tional mass increase with distance. Then the fundamental length scale i should be chosen 
in such a way as to ensure that the deviations from the Newtonian behavior begin at a 
certain distance from the gravitating body. A typical representative of such a situation is a 
spiral galaxy like our Milky Way, of mass Mg ~ lO^^M©, in which deviations from Newton's 
behavior (flat rotation curves) are prominent at a distance Vg ~ 10 kpc. This gives us the 
estimate 

£~ ~ 10 Mpc, (79) 
V rs 

where = 2GMg ~ 10~^pc is the Schwarzschild radius associated with the galaxy mass 
contained inside the radius Tg. Interestingly, this estimate for £ roughly corresponds to the 
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scale on which the relative perturbation Sp/p of the density in the universe becomes of order 
unity today. 

The critical distances Tc ~ {i'^Ts) for the Sun, for the Earth, and for an isolated proton 
are then given by 

r© ~ 2 pc , r© ~ 2 X 10'^ au , ~ 5 mm , (80) 

respectively. The large practical value of the critical radius ( 180|) even for microscopic bodies 
such as a proton will result in the effect that, in most (if not all) practical cases in the given 
example, the overall redshift factor (!63l) of low-curvature regions will be unobservable. In 
other words, the photons (and other particles) will always be emitted from regions of rela- 
tively high curvature (magnitudes of curvature corresponding to 3> 1 in the spherically 
symmetric case) and observed also in regions of high curvature, and there is no relative 
redshift factor between such regions. 

To explain this in more detail, assuming that we have gravitational sources separated 
by distances large compared to their critical scales ri*'' = ^f^ri^^j determined by their 
small-distance Schwarzschild radii ri^\ i = 1, . . . ,N, we can write the metric in the region 
far from all such sources as 

goo = i-J2 — , n:^rf, 2 = l,...,iV, (81) 

where rj is the distance from the point of observation to the z**^ source, and f't^ = Zri^^ 
is the renormalized Schwarzschild radius with the universal factor Z. We have used the 
superposition principle clearly valid in the domain i'^P <^ 1 and neglected the contribution 
from the effective cosmological constant. As we approach any of the sources within its 
critical radius, the metric is dominated by this source and reads 

gS = z'(^i-'-pj, n « r« , (82) 

with the same universal redshift factor Z. Because of this universality, a photon or any 
other particle traveling from one such source to another one will experience no additional 
relative redshift. 

Thus, for any sufficiently massive practical observer, only the effect of the apparent 
increase of the central mass (or, equivalently, of the gravitational constant) with distance 
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will be detectable, which asymptotically will result in the difference between and fg = Zr^. 
Although the notion of physical metric may not be well defined at intermediate distances 
(17^ . it is clear that the effect of apparent increase of the gravitational mass of the central 
body should be continuous with distance; hence, the value of the effective mass in our 
example will monotonically interpolate between its asymptotic magnitudes and Zr^ as 
one moves away from the gravitating body. This potentially can be used to explain the 
effect of missing gravitating mass in the universe. 



B. Pioneer anomaly 

It is interesting to see whether the theory under consideration might be able to explain 
the observed anomalous acceleration of the Pioneer spacecraft ap ^ 8 x 10^^'^ m/s^ [9,]. To 
explain the alleged missing mass in galaxies one has to take the function such that 
the effective gravitating mass increases with the distance. This would imply the presence of 
anomalous acceleration directed towards the Sun, potentially matching the observed value. 

To derive the exact observable result, we need to consider concrete realistic functions 
as well as the couplings of matter to the gravitational degrees of freedom of this theory, 
which will possibly distinguish one of the effective metrics from the whole conformal class. 
Some preliminary order-of-magnitude estimates can already be made irrespective of such 
details. The leading terms in non-metric corrections to the acceleration at small distances 
are expected to be of the order 

«o - , (83) 

as can be seen from the last terms on the right-hand side of flB4|) or flB5|) of Appendix [B] for 
the simplest generic case p2 = 1 [see also the concrete example (TfOl) ]. (To avoid confusion, 
we remind the reader that the speed of light was set to unity throughout this paper.) For 
the value of Ij^fa ^ 20 Mpc, the Pioneer acceleration ap ~ 8 x 10^^° m/s^ will be reached 
at a distance from the Sun r = 20 au. This value of ij^fa is of the same order as that 
which we obtained in the explanation of the missing mass in galaxies [see Eq. fl79|) ]. Thus, 
together with accounting for the missing mass in galaxies, the present theory may also be 
able to explain this observed anomaly in a uniform way. 

Because of the quadratic dependence on the distance in ( l83ll . the anomalous acceleration 
of the inner planets of the solar system will be well below the observational limits. However, 
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Object 


Distance (au) 


^thcor (10-10 m/s2) 


^obs (10-10 m/s^) 


Mercury 


0.39 


0.003 


0.04 


Icarus 


1.08 


0.02 


6.3 


Mars 


1.52 


0.04 


0.1 


Jupiter 


5.2 


0.5 


0.12 


Uranus 


19.2 


7 


0.08* 


Neptune 


30.1 


17 


0.13* 



TABLE I: The final column is the upper limit on constant acceleration determined from planetary 
orbits and taken from [l^. The constraints imposed by the orbits of Uranus and Neptune are 
somewhat uncertain; for this reason, they are marked with asterisk (see |lo|). The third column 
is the theoretical anomalous acceleration estimated from (j83p with ^/-y/a = 20 Mpc. Theoretical 
estimates for Uranus and Neptune, if proved to be valid, significantly exceed the observational 
constraints, which are also inconsistent with the Pioneer anomaly as a modification of gravity [lo| . 

the theoretical prediction for the outer planets based on the simple estimate ( !83|l seems to 
disagree with the observational constraints, as can be seen from Table [B The observational 
constraints for Uranus and Neptune, if proved to be valid, are also inconsistent with the 
general explanation of the Pioneer anomaly as a modification of gravity [lo| . 

The anomalous acceleration fl83l) is obtained for functions which have logarithmic 
asymptotic behavior at relatively large values of /3, as in example which results in 

metric fITOl) at small distances. If the function has power-law behavior at large values 
of /3, as in example (175]) . then the modified metric at small distances takes the form (1771) . 
leading to the anomalous acceleration of the order 

(84) 



a r 



ao - - 

r yi^Ts 

which differs from (1831) by a small factor r^/f^rg, of the order 10~^^ for Uranus and Neptune. 
In this case, the expected anomalous acceleration is many orders of magnitude smaller than 
that observed for the Pioneer, and satisfies well the solar-system constraints derived from 
planetary motion. 

Thus we can see that the specific prediction for anomalous accelerations within the solar 
system is strongly sensitive to the form of the unknown function 0(/3) in the domain of 
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relevant curvatures (3. It looks possible to explain the Pioneer anomaly (if it has gravitational 
origin) in frames of the modified theory of gravity under consideration, but a concrete 
realization of this requires more work. 



C. Non-cosmological redshifts of quasars and gamma-ray bursts 

The idea that high redshifts of quasars are not of cosmological nature, but rather are 
caused by some poorly understood physical circumstances, is pursued by a number of astro- 



physicists ll| (see also [12] for critical assessment and 13|] for a historical review). Similar 
considerations exist concerning the high redshifts of gamma-ray bursts One of the 

problems with this controversial hypothesis is to propose a viable alternative explanation 
for such high redshifts. Here we would like to demonstrate that our theory may be capable 
of providing such an explanation. 

Suppose that the condition > is satisfied in the range /?i < /5 < so that the 
redshift factor Z(/5i,/?2) in (I62!) is bigger than unity. Then, assuming for the moment that 
the end points of this range are in the metricity regime 0/? <^ 1, and taking into account 
that the value of j3 in the corresponding regions is the value of the Weyl curvature of the 
physical metric, we conclude that light emitted from regions of higher Weyl curvature will be 
additionally redshifted by the factor Z(/5i,/32) > 1- The observable redshift ;Zobs of a source 
with such intrinsic redshift factor Z and with local (peculiar or cosmological) redshift z will 
then be given by 

^obs = Z{l + z)-l. (85) 

This effect can be responsible for the observed high redshifts of quasars and gamma-ray 
bursts. 

The value of Z{(3i,(32) is quite arbitrary depending exponentially on the behavior of the 
function 0(/9). As a simple example, consider the function 

0(/5) =0ilog(l + x2) , x = e(3, (86) 

which gives the "metric" behavior in the two domains x -C 1 and x ^ 1. The derivative 
(f)f3 should also satisfy the condition < 2 in order that we do not encounter a singularity 
outside the horizon. The validity of this condition for all (3 leads to the constraint 

-=^<l. (87) 
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The redshift factor is given by the expression 

a;2>l 




Z = exp I / I ^ e™ < 23.14 , (88) 



2x 

where the last upper bound comes from the bound in ( IHTI) . The upper bound ( IHTI) . ( IHHjl 
for the admissible values with function fl86l) is more than enough to explain the redshifts of 
quasars and gamma-ray bursts up to ^obs ^ 6. Indeed, we need a ~ 0.62 to get the highest 
values of Zobs ~ 6 for 2; = 0. 

According to the described scenario, additional redshift factor Z will be present in all 
massive compact objects in which radiation originates in regions with curvature (3 > 
i.e., having the following relation between the size R and Schwarzschild radius : 

R'<rsi'. (89) 

The free parameters i and a should be fitted to model redshifts of compact objects with 
strong gravity. Thus, if we take an accretion onto a black hole as a working model for a 
quasar, then we will have ~ r^, and condition fl89l) will become Vg < i, implying that only 
black holes of mass smaller than some value produce redshifted emission. To fix the value of 
i, one needs to build a detailed model of a quasar, which obviously goes beyond the scope of 
this paper. We only stress that now it has to be built under the assumption that quasars are 
not superluminous objects situated at very large distances but rather are compact massive 
objects of average absolute luminosity in our proximity. Also note that the gravitational 
redshift in our theory is different by nature to the usual gravitational redshift from a massive 
body in GR. This is clear from the observation that our effect involves a length scale £ in 
( I89ll and, therefore, can take place at arbitrary Weyl curvatures, for example, at which the 
usual redshift would be negligible. The redshift itself depends on another free constant a, as 
can be seen from f lHHj) . Moreover, our effect in a wide range of curvature values around the 
central body can be almost independent of the value of curvature itself. Hence, the usual 
constraints jl5(] that rule out the gravitational redshift of quasars in general relativity have 
to be revised in the theory under investigation. 

There are two possible sources of the observed scatter in the values of redshifts. The 
redshift factor Z{Pi, P2) will depend on the actual curvature j32 of the region where radiation 
is formed. We remember that the effective physical metric is not yet defined in the region 
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where |0/3| ~ 1; however, the observable effect of redshift from these regions is expected to 
be of the magnitude Z{Pi,P2)- The spread in the values of will then lead to a spread in 
the observed redshifts. Another obvious source of scatter in the redshift values is the overall 
cosmological expansion and peculiar motion which contribute to the observable redshift 
(l85ll . As is clear from (l85l) . the effect due to peculiar motion and cosmological expansion is 
amplified by the factor Z. For example, if Z = 2, then a nearby quasar situated at z = 
will be observed at Zobs = 1, while a quasar situated at 2; = 0.5 will be seen to have Zobs = 2. 

The non-cosmological redshift effect which we consider in this subsection could also be 
relevant to neutron stars. In order that this effect be appreciable, condition ( l89l) should 
be satisfied at their surfaces. Considering that neutron stars have approximately nuclear 
density, we conclude that the ratio R^/vg (where R is their radius) is approximately constant 
in such objects, and condition fl89l) then gives a numerical estimate for i : 

£ > 30 km . (90) 

With i satisfying this condition, radiation from the surfaces of neutron stars will be addi- 
tionally effectively redshifted. This effect, in principle, may be tested by observations. 

Concluding this subsection, we note that the same condition f l90|) will formally imply 
nontrivial gravitational effects from the region of atomic nuclei. However, the classical 
picture of gravity that we considered here may be not valid in such microscopic regions. 
This question deserves additional investigation. 

D. Combined scenario 

In order to implement simultaneously the described physical effects, the function 0(/3) 
much be chosen appropriately. It is clear that we require the existence of three different 
domains in the (3 space, with solar system in the region characterized by the "metricity" 
property 0/3^1. Let us denote representative values of P in the order of increasing by Pi, 
(32, and jS^. The solar-system values of curvature will correspond to (32- We should have 
Z{(3i,f32) < 1 to explain galactic rotation curves, and Z{j32,(33) > 1 to explain high redshifts 
of quasars and gamma-ray bursts. The function that may do the job will look something 
like 

(j)i(3) = (/.o - log (1 + Xl) + 02 log (1 + Xl) , Xi = il(3 , X2 = il(3 , (91) 
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FIG. 2: Typical qualitative shape of the function (f){f3) which might exhibit the three effects dis- 
cussed in this paper, namely, missing mass in galaxies, anomalous acceleration in the outskirts of 
the solar system, and high redshifts of quasars and gamma-ray bursts. The titles of astrophysi- 
cal objects are placed at the corresponding values of /?, which is roughly the Weyl curvature of 
space-time characteristic of these objects. 

with ii ^ £2 and <^2 > The value of <^o should be negative to account for the positive 
large-scale cosmological constant A = — 30o- Qualitative behavior of the function is 
shown in Fig. O However, one can propose other functions with the desired properties, 
and it is by combined theoretical investigations in quantum gravity and by comparison with 
observations that a reasonable function (f){(3) may be established. 

VII. SUMMARY 

In this paper, we considered the vacuum spherically sjTumetric problem in the theory 
of gravity described by action ([1]), which modifies the Plebahski self-dual formulation of 
general relativity by introducing an additional function (p [Tr (\I'^) , (^I/^)]. In the spherically 
symmetric case, the two arguments of the function are not independent, and it becomes 
a function of only one effective field (3, defined in ([6]). The form of the function is to be 
determined from quantum considerations. Since this problem is not yet solved, in concrete 
physical applications we have explored several different possibilities of its form. 

The basic variables of the theory under consideration are the spin connection one-form 
A and the two-form spinor field B, while the Lagrange multiplier field can be thought 
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of as expressible in terms of B via the algebraic "metricity" equations. The theory with 
nontrivial does not distinguish any special metric; however, there arises a conformal 
class of metrics with respect to which the two-form B is self-dual. The existence of such a 
conformal class enables us to speak about the conformal structure of the solution. 

The vacuum spherically symmetric solution has several significant features, and we can 
summarize our results as follows: 

1. For arbitrary (p, the theory respects the analog of the Birkhoff theorem, i.e., the vacuum 
spherically symmetric solution is necessarily static and asymptotically describes a space of 
constat curvature. In our theory, this means the existence of a vector field generating a group 
of diffeomorphisms that leave invariant the fields A, B and In the coordinate language, 
the components of these fields do not depend on the "time" coordinate t. This property 
is probably connected with the fact that our modifications of the Plebahski formulation of 
the general-relativistic action does not increase the order of differential equations. A similar 
property was recently proved for metric theories which preserve the second-order character 
of field equations 

2. In any domain of its argument in which the function (f){(3) varies slowly (0/3 <^ 1), 
it acts simply as a (multiple of) the cosmo logical constant, and the solution in the corre- 
sponding spatial domain possesses approximate "metricity" property in the sense that the 
three components of the self-dual spinor two-form B are expressible as self-dual parts of 
the canonical exterior products of some basis one-forms. These basis one-forms then define 
a unique metric in the corresponding space-time region, which is an approximate solution 
of the vacuum Einstein equations. In our spherically symmetric solution, such a behavior 
is obtained at spatial infinity, where we recover the Schwarzschild metric with small and 
rapidly decaying "nonmetric" corrections. The value of j3 in this case plays the role of Weyl 
curvature. 

3. If several regions exist in the domain of f3 where the function (f){(3) is slowly vary- 
ing, then the distinguished metric in the corresponding spatial regions is described by the 
Schwarzschild-de Sitter form but with different values of the Schwarzschild radius and ef- 
fective cosmo logical constant. In other words, the gravitational and cosmological constants 
become curvature-dependent in this theory. This property of the solution can potentially be 
used to account for the problem of missing mass in spiral galaxies and other astrophysical 
objects. The nonmetric corrections in the regions close to the gravitating bodies in this case 
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could also explain the observed anomaly in the acceleration of the Pioneer spacecraft. 

4. In addition to this effective remormalization of gravitational and cosmological con- 
stants, there arises a nontrivial universal redshift factor between regions of different Weyl 
curvature (3. This effect potentially can explain the high redshift of quasars and gamma-ray 
bursts. 

5. The conformal structure of our solution inside the black- hole region is different from 
that of the Schwarzschild solution and depends on the form of the function cj). Typically, 
"inside" the analog of the Schwarzschild horizon one finds another surface of extreme "non- 
metricity" . This surface is spacelike, and replaces the usual spacelike singularity inside the 
Schwarzschild black hole. The metric ceases to exist at this surface, but all the dynamical 
fields of the theory are finite. Thus, this surface is only a metric singularity, but not a 
singularity of the theory. Across this surface, the coordinates t, r change their spacetime 
roles once more, and one typically finds another Schwarzschild-like horizon behind this 
"non-metricity" surface. The theory does not admit arbitrarily small black holes: for small 
objects, the curvature on the would-be horizon is so strong that non-metric modifications 
prevent the horizon from being formed. Instead of horizon, one has "naked" hypersurface 
of non-metricity in this case. The details of the conformal diagram depend on the specific 
shape of the function For a simple choice = it is shown in Fig. [H 
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APPENDIX A: PROOF OF THE STATIC PROPERTY OF THE METRIC 

In this appendix, we give the details of the proof of static property of the solution to the 
vacuum spherically symmetric problem. 

The most general spherically symmetric expression for the two-form B is given by (HM : 



B = X_c m Al + X+c nAm + X{lAn — mA m) 



(Al) 



with the one-forms /, n, m, and m given by (|T6 



/ = — {fdt - gdr) 



1 _ r 

n = — ^ {fdt + gdr) , m, ffi = —j= [hdO ± i sin 6 d(j)) . (A2) 



Now the functions /, h, and c are not assumed to be time- independent, but are functions 
of both r and t. We are going to show that h is a constant (coordinate-independent), and 
f,g, and c are time-independent due to the field equations. 

The expressions for B are now modified because of the appearance of time derivatives, 
so that instead of fl20l) . we have 



dB. 



{rcfY 



V2r 

, (eg)- 
V2fg 



gf 



I An Am — c cot 91 Am Am 



TlTTlh 

I An Am + - — ^ _ I An A {m + m) 



dB, 



gf 



dBn 



V2r 

I jcgy 

V2fg 

[I — n) A m A m 



2V2fh 

I A n A m + c cot 9 n Am Afh 



(A3) 



I An Am + 



ch 



2V2 fh 
h 



I An A {m + m) , 
{I + n) Am Afh . 



rg ' ■ V2fh 
where overdot denotes the time derivative. Formulas fl22l) are then modified as follows: 



(dS_ A m) 
{dB+ A m) 



ch 



2V2 fh 
{rcfY (eg) 



ch 



V2rgf V2fg 2^2 fh 



, {dB+ A m)^ 



ch 



2V2 fh 



(A4) 



- (rf5_ A n)^ = {dB+ A /), 



V2'i 



cot 9 . 



{dB^An)^ = ^ 



h 



rg V2fh 



{dBo A /),„ = — + ^ 



rg V2fh 
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the rest of the projections being zero. Then, using formulas (125|) and (1261) . which remain 
unchanged, we obtain 

A_ = c^^Pm + cQfh , = cQm + c^^Prh , 

cot 6* h (cgVh) ^^^^ 

where we made the notation 

1 (r^h)' h 



P 



2V2g r^h 2V2fh 



h' h 



R 



2V2gh 2V2fh 
rcfVh 



\f2rcf g\fh 

Now we have to compute the left-hand sides of equations fl3Tl) and equate them to zero. 
Computing the m A m component of any of the first two equations in ( |3Tll . we immediately 
obtain that Q = 0, which, in turn, simplifies the expressions ( lASp . Calculating then the 
(/ A fn) and n Nfn components of the first equation in fl3T]) . we get 

respectively. The condition P = is excluded because it contradicts the / A m component 
of the same equation, which contains a generically nonzero expression (3 + (j). Therefore, one 
must have h = const. 

The remaining four equations stemming from the first two equations in ( 13T|) read 



V2rg V2cfg, ^ v^/ 



V2rg V2cfg, ' ' V2 f 

+ ^ + c-P (c-P + P) + ^ -ciP + 0) 



(A8) 



V2rg V2cfg. ' ' V2 f 



V2rg V2cfg, ^2 f 

where = eg, as usual. Now, subtracting the first equation from the third one, and the 

second from the fourth one, we obtain, respectively, 

:^ + (c-ip)=o, :^-(c-ip)=o, (A9) 

eg* ^ ' eg* ^ ' 
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which imphes 

()* = 0, (c-ip)=0. (AlO) 

These two equations are equivalent in view of the condition h = const and definition ( ]A6p . 

Under condition (lAlOl) . the remaining two equations stemming from (1A8I1 are precisely 
equations ( l32l) . while the last equation in (l3Ti) leads to the two equations ( l33l) . Thus, we 
have the system of differential equations fl35|) . The second equation in fl35l) then implies that 
P does not depend on time, hence, by virtue of (fT5l) . c is also time-independent. Finally, 
the third equation in fl35|) implies that / can only have a time- dependent overall factor, 
which can always be rescaled to a constant by changing the time variable. This completes 
the proof of the static property of the spherically symmetric vacuum solution in the theory 
under investigation. 



APPENDIX B: MASS RENORMALIZATION AND REDSHIFT: GENERAL 

ANALYSIS 



In this appendix, we generalize the examples considered in Sec. |Vl Thus, assume that 
two regions are characterized by the conditions £"^(3 <^ 1 and i'^jS ^ 1, respectively, where i 
is some length scale. We will find the approximate solutions in these regions assuming the 
behavior 

2aixP^ + o(xPi) , a; < 1 , 
2a2X-P^ +o{x-P^) , a; > 1, 

where the variable x = £"^13, and ai, a2, pi > and ^2 > are different constants. 
At small radial distances, where x ^ 1, we integrate equations fj42l) to obtain 



2r3 



X 



1 + —x-P' + o (x~P') 

P2 



f^g^ = 1 + ^x-P' + O (X-P^) , 

X > 1 



(Bl) 



(B2) 



We also have 

c2 = 1 - 3a2x'P^ + o {x~P^) , X > 1 . (B3) 
Collecting all the terms, we obtain the leading contribution to the goo coefficient in metric 

mi- 

(B4) 



^700 = /' = 1 - ^ + 0oor' + 3«2 ( 1 + 77- 



3P2; V^rJ 
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where = const is the value of the effective Schwarzschild radius at small distances, 0oo 
is the value of the function at infinity. We have also assumed the condition <^ r, 
which guarantees that some terms that are of lower power in r are actually subleading. 
This expression is valid if p2 7^ 1- In the interesting case p2 = 1, the correction to the 
Schwarzschild-de Sitter metric will also have logarithmic terms; see our concrete example 
(IMD in Sec. |Vl 

As we discussed above, the physical metric can be specified by an additional conformal 
factor, in which case, we obtain somewhat different corrections. Thus, for the Urbantke 
metric defined in ( l52l) and for the "volume" metric defined in (Ell), we will have the coinciding 
approximate relations for -C i'^Vs : 



+ + 2^2 1 + 



goo = i 

At large radial distances, where x <^ 1, we have 

P2^ 



P2 



2^3 



P2 



(B5) 



exp 



2r3 
f*g* = exp 







2/3 



dp 



Pi 



1- —xP' + o{xP'' 

Pi 



c2 = 1 - Saix""' + o {x^') , X < 1 . 
In this case, we obtain the following leading contribution to the metric: 



where 



1 - 



0(0), and 



+ 0or + 3ai 1 



3pi 



2r3 



Z = Z(oo,0) = exp(^-^ ^dp^ 



(B6) 
(B7) 

> e^Zvs , (B8) 
(B9) 



In the case of extra conformal factors, again, the corrections to the Schwarzschild-de Sitter 
metric will be somewhat different. For the Urbantke metric defined in ( |52l) and for the 
"volume" metric defined in flMl) , we have the same approximate relations for ^ ^^r^ : 



fi'oo 



1 - 



+ <t)y + 2ai 1 

Pi 



l^Zr, 
2r3 



(BIO) 



We note that ^ 1 in our approximation in these regions, so the physical metric is 
well defined. The contributions proportional to 02 and ai in (1B4I) . (IBSP and (1B8I) . (1B10I) . 
respectively, are small and can be dropped if the constants ai and 02 are not very large. 
However, these contributions themselves are of the same order as the "nonmetricity." Their 
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physical interpretation, therefore, requires the knowledge of the matter couplings in our 
theory, which is an issue still to be resolved. 

For Z > 1, the observed gravitational mass of the central object at large distances is Z 
times larger than it is at small distances. In the intermediate region of radial distances 

1<^^^> (Bll) 

the metric coefficient behaves in a complicated way; moreover, the physical metric may not 
be well defined in this region at all. If several regions exist in the space of (3 in which <^ 1, 
then obvious renormalizations of the observed gravitational masses and redshift factors exist 
between these regions. 
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